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O ' Abstract 
D ■ 

Q , We apply Bayesian methods to the study of the sensitivity to neutrino masses of a 

I Galactic supernova neutrino signal. Our procedure makes use of the full statistics of 

events and is remarkably independent of astrophysical assumptions. Present detectors 
can reach a sensitivity down to rriu ~ 1 eV. Future megaton detectors can yield up 
I to a factor of two improvement; however, they will not be competitive with the next 

O ■ °' ^-'^^^^ '^-'^^^^ experiments. 
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^ ' 1 Introduction 
p , 

It was realized long time ago that neutrinos from a Supernova (SN) can provide valuable 

Qh, informations on the neutrino masses. The basic idea relies on the time-of-flight (tof) 

Cl^l delay At that a neutrino of mass niiy and energy E^^ traveling a distance L would suffer 

^ ' with respect to a massless particle: 

^ r / L \ /10MeV\y m^ \2 
At = 5.1ms — ]\—^ ^7 • (1) 
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>y , Indeed the detection of about two dozens of neutrinos from SN1987 allowed to set 

model independent upper limits at the level of m-p^ < 30 eV Since SN1987A, sev- 
eral proposal have been put forth to identify the best ways to measure the neutrino 
tof delays. Often, these approaches rely on the identification of "timing" events that 
are used as benchmarks for measuring the time lags, as for example the simultaneous 
emission of gravitational waves or the abrupt interruption of the neutrino flux due 
to a further collapse into a black hole [S]. The less model dependent limits achievable 
with these methods are at the level of rrii, ^ 3eV, and tighter limits are obtained 
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only under specific assumptions. A different method to extract informations on the 
neutrino masses that uses the fuh statistics of the signal and does not rely on any 
benchmark event was proposed in Ref. 4 and developed in Ref. [3]. Here we resume 
the results obtained by applying this method to determine the sensitivity of the Su- 
perKamiokande (SK) water Cerenkov and KAMLAND scintillator detectors, and of 
the planned experimental facilities Hyper-Kamiokande and LENA. 



In real time detectors, supernova i^g are revealed through to the positrons they produce 
via charged current interactions, that provides good energy informations as well. Each 
Ve event corresponds to a pair of energy and time measurements {Ei,ti). In order to 
extract the maximum of information from a high statistics SN neutrino signal, all the 
neutrino events have to be used in constructing a suitable statistical distribution, as 
for example the Likelihood, that can be schematically written as: 



Ci represents the contribution to the Likelihood of a single event with the index i run- 
ning over the entire set of events, (t{E) is the detection cross-section and F{E; t) is 
the energy spectrum of the neutrinos, whose time profile can be reconstructed rather 
accurately directly from the data [H [3] . The main problem in constructing the Like- 
lihood ((21) is represented by the (unknown) time profile of the neutrino flux (j){t). We 
construct a flux model by requiring that it satisfies some physical requirements and 
a criterium of simplicity: i) the analytical flux function must go to zero at the origin 
and at infinity; ii) it must contain at least two time scales for the neutrino emission 
corresponding to the fast rising initial phase of shock-wave breakout and accretion, 
and the later Kelvin-Helmholtz cooling phase; iiij it must contain the minimum pos- 
sible number of free parameters. The following model for the flux has all the required 
behaviors: 



The five parameters that on the l.h.s of Q have been collectively denoted with A are: 
two time scales ta for the initial exponentially fast rising phase and tc for the cooling 
phase, two exponents Ua and ric that control their specific rates and one additional 
exponent Up that mainly determines the width of the "plateau" between the two phases. 
Since we are interested only in the neutrino mass squared m^, irrespectively of the 
particular values of the nuisance parameters A, starting from the Likelihood ^ we will 
need to evaluate the marginal posterior probability p{m^\D), that is the probability 
distribution for mf, given the data D. This is done by marginalizing the posterior 
probability with respect to the nuisance parameters: 



2 The method 
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p{m1, A|/), that in Bayesian language is called the prior probability of the model, allows 
us to take into account any available prior information on the parameters and A. 
We will use flat priors for all the A's and, to exclude unphysical values of m^, a step 
function @{ml) = 1, (0) for ml > 0, (< 0). 

The dependence on could be directly included in the flux ^ by redefining 
the time variable according to However, it is more convenient to proceed as 

follows: given a test value for the neutrino mass, first the arrival time of each neutrino 
is shifted according to its time delay, and then the Likelihood is computed for the 
whole time-shifted sample. Subtilities in the evaluation of the Likelihood contribution 
Ci{ti,Ei) arising from the uncertainty AEi in the energy measurement are discussed 
in Refs. gUg. 

We have tested the method applying it to a large set of synthetic Monte Carlo 
(MC) neutrino signals, generated according to two different SN models: SN model I 
corresponds to the simulation of the core collapse of a 20 Mq star ^6 carried out by 
using the Livermore Group code . In this simulation I'^^r opacities were treated in a 
simplified way, and this resulted in quite large (and probably unrealistic) differences in 
their average energies with respect to i>e . SN model II corresponds to a recent hydro- 
dynamic simulation of a 15 Mq progenitor star |7^ carried out with the Garching group 
code jU] • This simulation includes a more complete treatment of neutrino opacities and 
results in a quite different picture, since the antineutrino spectra do not differ for more 
than about 20%. In both cases the effects of neutrino oscillations in the SN mantle 
have been properly included in the simulations. Note that the two types of neutrino 
spectra of SN model I and II fall close to the two extremes of the allowed range of 
possibilities. This gives us confidence that the results of the method are robust with 
respect to variations in the spectral characteristics. 

3 Results 

To test the sensitivity of our method, we have analyzed a large number of neutrino 
samples, grouped into different ensembles of about 40 samples each. For each ensemble 
we vary in turn the SN model (model I and II), the SN-earth distance (5, 10, and 15 kpc) 
and the detection parameters specific for two operative (SK and KamLAND) and two 
proposed (Hyper-Kamiokande and LENA) detectors. The results of our recent detailed 
analysis [5^ are summarized in Tabled In columns 2 and 4 we give the 90% c.l. upper 
limits that could be put on mj, in case its value is too small to produce any observable 
delay. In columns 3 and 5 we estimate for which value of rrii, the massless neutrino case 
can be rejected at least in 50% of the cases. These results confirm the claim ^ that 
detectors presently in operation can reach a sensitivity of about 1 eV, that is seizable 
better than present results from tritium /3-decay experiments, competitive with the 
most conservative limits from neutrinoless double /3-decay, less precise but remarkably 
less dependent from prior assumptions than cosmological measurements. However, in 
spite of a sizeable improvement, future detectors will not be competitive with the next 
generation of tritium /3-decay and neutrinoless double /? decay experiments. 
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1 9 


b) SK (5 kpc) 






1.1 ±0.3 


1.0 


c) SK (15 kpc) 








1.6 ±0.6 


1.4 


d) SK+KL (10 kpc) 


1.0 ±0.2 


0.9 


1.1 ±0.3 


1.0 


e) HK (10 kpc) 


0.4 ±0.1 


0.4 


0.5 ±0.1 


0.5 


f) LENA (10 kpc) 


0.9 ±0.2 


0.9 


0.9 ±0.3 


0.9 



Table 1: Results for the fits to m^: a)-c) SK for different SN distances, d) SK plus Kam- 
LAND, e) Hyper-Kamiokande, f) LENA. All masses are in eV. 
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